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Soot formed in laminar diffusion flames of heavily sooting fuels evolves through four distinct growth 
phases which give rise to four distinct aggregate fractal morphologies. These results were inferred from 
large and small angle static light scattering from the flames, microphotography of the flames, and analysis 
of soot sampled from the flames. The growth stages occur approximately over five successive orders of 
magnitude in aggregate size. Computer simulations can reproduce the kinetics observed in aggregation 
and aerogelation of soot clusters and indicate that these four growth stages involve either diffusion limit 
aggregation or percolation in either three or two dimensions. 
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Abstract 

Soot formed in laminar diffusion flames of heavily sooting fuels 
evolves through four distinct growth stages which give rise to four 
distinct aggregate fractal morphologies. These results were inferred 
from large and small angle static light scattering from the flames, 
microphotography of the flames, and analysis of soot sampled from the 
flames. The growth stages occur approximately over four successive 
orders of magnitude in aggregate size. Comparison to computer 
simulations suggests that these four growth stages involve either 
diffusion limited cluster aggregation or percolation in either three or two 
dimensions. 
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Fig.l 


Diagram of the small angle and large angle 
light scattering apparatus. These apparatii 
allow us to probe soot morphologies from 70 
nm to nearly 60 pm. 
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Fig. 2 Scattered Intensity as a function of wavevector, I(q), for 

an acetylene flame for various heights h above the burner. 


Note the presence of: D — 1.8 aggregates with size ~ 0.3 JLL m and 

D —2.6 aggregates with size ~ 12 JLL m at large h. 
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Right side of Fig.2b shows an example of submicrom soot (Stage 1) 
while left side shows supermicron soot superaggregates (Stage 2). 

Similar behavior was observed for Toluene and 1-methylnapthalene. 
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Fig. 3 Transmission electron micrographs of large soot 
(supersoot) aggregates. 

Next to each image is the perimeter of the aggregate. 

Values for the radius of gyration R g and mass fractal 
dimension D are given for each cluster. D was 
determined from the perimeter fractal dimension D p 
following Julien, Thouy, and Ehrburger-Dolle (1994). 


TEM Analysis (D>2) 



R g =17juun Rg= 1 7.3[im R g -14.4^m 

D - 2.53 D ~ 2.51 D-2.48 
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Fig. 4 Structure factor of an off-lattice, DLCA simulation in a cubic 
box (3d) of side L (measured in primary particle radii a) with 
N m primary particles. 

The structure factor S(q) is proportional to the light 
scattered intensity I(q) for a real system. 

S(q) evolves with increasing time, as represented by the 
decreasing number of clusters, N c , from a system with solely 
3d, DLCA, D - 1.8 aggregates to a system with percolated 
superaggregates with D - 2.6 composed of the smaller 
D =d.8 aggregates. 
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Fig. 5 Diagram of a Superaggregate 


Superaggregate 

An aggregate of one morphology composed of 
smaller aggregates of a different morphology. 



length scale, 
(superaggregate). 
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Fig. 6 (a) 7 nsec exposure photograph of soot in an 
acetylene/air flame at a height above burner of 
h=10cm. The width of the picture is 1mm. 

(b) Binary digitized version of (a). 

(c) Number of pixels per aggregate versus radius 
of gyration for an ensemble of soot clusters 
including those in (a). 
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Fig. 7 Fourier transform, hence structure factor S(q), of 
the digitized image of a flame soot network, lower 
left corner, in a laminar acetylene/air diffusion 
flame at a height above burner of h=12cm. 

Note the crossover from a slope of - 1.4 to a 
slope of - 2.1 indicating the formation of a 2d 
aerogel network. 
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Fig. 8 Structure factor of an off-lattice, DLCA simulation on a 
square lattice (2d). 

The q' 3 0 regime is the Porod regime of the square box with 
q uMl , d=2 and interference ripples due to the sharp edges of 
the box. 

The q' 21 regime is due to percolated superaggregates with a 
fractal dimension of 1.9. 


The q" 1 ’ 45 regime is due to 2d, DLCA aggregates that make up 
the superaggregate. The volume fraction of monomers is 
f v =0.25. 
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CONCLUSIONS 


• Current work leads to new picture of Soot 
Morphology in Heavily Sooting Laminar 
Diffusion Flames. 


• With increasing Soot aggregate size, a rich 
phenomenology unfolds. 

• Four different stages of growth — each 
causing a specific morphology in the 
resulting aggregates. 

• Summary of Current Work: 


Table I. The Realm of Supersoot 


Stage 1 
3d DLCA 


Stage 2 Stage 3 

3d Percolation or 2d DLCA 
Restructuring? 


Stage 4 

2d Percolation 


£> = 1.8 D-2.6 £> = 1.4 £> = 1.9 

aggregates superaggregates supersuperaggregates gel network 




~ 10(1/77 


~ lOOjim 
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The current view of soot formation and growth in flames involves a series of steps including fuel thermal 
decomposition to small radicals (Palmer and Cullis , 1965; Glassman, 1988) that react to form polyaromatic 
hydrocarbons (D’Anna, D’ Ales so and Minutolo, 1994; Dobbins and Subramaniasivam, 1994; Dobbins, Fletcher 
and Chang, 1998; Frenklach, 2002) that subsequently nucleate, coalesce and dehydrogenate to yield roughly 
spherical, graphitic , primary particles of soot with diameters in the few tens of nanometer range (Lahaye and Prado, 
1981). After this, the physical process of three-dimensional diffusion limited cluster aggregation (DLCA) (Meakin, 
1999; Oh and Sorensen, 1997) proceeds to make noncoalesced clusters which have a fractal morphology with a 
universal fractal dimension of D ~ 1.8 (Samson, Mulholland and Gentry, 1987; Zhang, Sorensen, Ramer, Olivier 
and Merklin, 1988; Dobbins and Megaridis , 1990; Sorensen, Cai and Lu, 1992; Koylu and Faeth, 1992). Thereafter 
environmental restructuring might occur, but with this single step of DLCA, the present description of in-flame soot 
growth ends. 

Our work of the past few years has indicated that this picture is seriously incomplete for the soot produced 
in heavily sooting laminar diffusion flames. This includes the observation of gel-like networks of soot (Sorensen, 
Hageman, Rush, Huang and Oh, 1998), flat, two-dimensional soot (Sorensen and Hageman, 2001), and soot with a 
fractal dimension of 2.6 that was found in a variety of different flames (Sorensen, Kim, Fry, Shi and Chakrabarti, 
2003; Kim, Sorensen and Chakrabarti, 2004). In this report we show that with increasing soot aggregate size, a rich 
phenomenology for soot in these flames unfolds. Four different stages of growth occur each causing a specific 
morphology in the resulting aggregates. These results are displayed in Table I (of main poster, conclusions section), 
and the evidence for these results will be described below. Our work is relevant to soot formation in common 
laminar diffusion flames in which the fuel and the oxidizer (e.g., oxygen in air) flow in a parallel manner in adjacent 
regions and diffuse together at the flame front. 

In our experiments static light scattering was used to measure the structure of the soot aggregates in the 
flames (Sorensen, 2001). These measurements gave the scattered intensity I(q) as a function of q, the scattering 

wave vector, q = 4nX~ l sin(0/2) where X is the optical wavelength and 0 is the scattering angle. An argon ion 
laser operating at X = 488nm was used as the light source. Both a small angle apparatus (Ferri, 1997), with 

150 cm" 1 < q < 2.7 x lO 4 cm" 1 (ca. 0.07° < 0 < 12° ) and a large angle apparatus, with 10 4 cm"' < q < 1.4 x 10 5 cm * 

(ca. 5 < 0 < 70° ) were used. A diagram of the apparatus is presented in Fig. 1 (of main poster). Since q 1 is the 
length scale of the scattering experiment, these apparatii allow us to probe soot morphologies from 70nm to nearly 
60 pm . 

Direct observation of the very large soot in the flame and the flame structure was accomplished using a 10 
power photomicroscope with an object distance of 15 cm. The flame was backlit by a 7ns pulsed Nd:YAG laser 
with X = 532nm shining on a white background behind the flame. The flames were all laminar diffusion flames in 
ambient air. Gaseous fuels at controllable flow rates emanated from a brass tube 0.9 cm inside diameter. Liquid 
fuels were burned with a simple wick burner with wick diameter of 0.6 cm. Solid naphthalene was molded into a 
candle 0.6 cm in diameter. The fuels used are listed in Table 2 below. 

Figure 2a of main poster shows the scattered intensity as a function of the wavevector, I(q), for an 
acetylene flame. Low in the flame, i.e., for heights above burner of h<2cm, the scattering is isotropic for low q but 

then bends over to a power law with slope ca. -1.8. The bend occurs near q ^ 3 xlO 4 cm" 1 to imply an aggregate 

size of q" 1 — 0.3)Lim . This corresponds to the often observed and well known phase of submicron, fractal dimension 
D — 1.8 soot which is a result of a 3d DLCA aggregation process (Meakin, 1999; Oh and Sorensen, 1997) (stage 1, 
Table 1). 
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Table 2. Fuels, threshold soot indices (TSI), and measured fractal dimensions within the 4 stages of soot 
growth. Dash means no measurements were attempted although stage 1 is expected for all flames. 


Fuel 

TSI 

Stage 1 

Stage 2 

Stage 3 

Stage 4 

Gases 






Methane 

~0 

1.8 

No 

No 

No 

Ethylene 

0.7 

1.8 

No 

No 

No 

Acetylene 

11 

1.7 

2.7 

1.44 

1.9 

Fiquids 






Hexadecane 

5.8 

- 

No 

No 

No 

Isooctane 

6.4 

- 

No 

No 

No 

1-Hexadecene 

11 

- 

No 

No 

No 

JP5 Fuel 

- 

- 

Marginal 

No 

No 

Decaline 

15 

- 

Marginal 

No 

No 

Toluene 

44 

1.7 

2.6 

1.39 

1.9 

Styrene 

67 

- 

2.6 

- 

- 

1 -Methyl 
Naphthalene 

91 

1.8 

2.5 

(1.4) 

1.9 

Solid 






Naphthalene 

100 

- 

2.7 

1.5 

1.9 


With increasing height above burner, which corresponds to increasing time for the soot to grow, a very 
intense, second phase appears at low q. In the region ca. 800 < q < 8000 cm 1 the slope of the log-log plot indicates a 

power law I(q) ~ q" 26 for this second phase. Note that the q 4 ' 8 regime remains at large q. This scattering behavior 
may be interpreted to indicate a phase of ca. 12pm soot with a fractal dimension ofD ~ 2.6 for length scales of 


1.2pm to 12pm (stage 2) and D— 1.8 for scales smaller than 0.3pm down to the primary particle size (ca. 25nmas 


determined by TEM ). Results for the toluene and 1-methylnaphthalene show the same behavior. We have shown 
(Kim, Sorensen and Chakrabarti, 2004) that this large, D — 2.6 phase occurs in laminar diffusion flames of all 
heavily sooting fuels as quantified by the Thermal Sooting Index (TSI) as defined by Calcote and Manos, 1983 and 
Olson, Pickens and Gill, 1985, specifically when TSI > 10 ± 2 (roughly, the larger the fuel TSI, the sootier the 
flame). Otherwise, for less sooty flames only the submicron, D — 1.8 soot is seen; no D — 2.6 phase nor the 
subsequent phases to be described below (see Table 2). 

To corroborate the light scattering measurements soot was sampled from the flame using a well established 
thermophoretic technique (Dobbins and Megaridis , 1987). Figure 2b (right side) shows an example of submicron 
soot (stage 1) collected low in the acetylene/air diffusion flame and Fig. 2b (left side) shows an example of 
supermicron soot superaggregates (stage 2) collected high in the flame . Visual inspection shows the much more 
compact appearance of the large soot than the small, qualitatively supporting a larger fractal dimension for the larger 
aggregate. Previous real space analysis of projected images of submicron soot have confirmed the D — 1.7 to 1.8 
fractal dimension (Dobbins and Megaridis, 1987; Koylu and Faeth, 1992; Cai, Fu and Sorensen, 1993). 

For the large soot aggregates we applied an analysis developed by Jullien, Thuoy and Ehrburger-Dolle, 
1996 for 2d projected images of 3d fractal aggregates with D>2. Using simulations they found the relation 


/ \ 2 / 3 

D = 3- ( D p -1J between the perimeter fractal dimension D p and the true 3d mass fractal dimension D. Figure 3 


in the main poster gives more examples of these large, relatively dense soot aggregates. We measured 10 large 


clusters to find D p = 1.26 ± 0.06 . This D p implies D = 2.6 ± 0.06 , in excellent agreement with our light scattering 


results. Furthermore, the average size we measure for the collected aggregates is 10 ± 3pm , also in excellent 
agreement with light scattering. 

The mechanism that creates this second phase is still somewhat uncertain. One scenario is that it represents 
3d DEC A clusters with initial D — 1.8 that have grown to ca. lOjLtm or larger and then restructured to a denser 

D ~ 2.6 (hence smaller) under the shear stress at the flame front. There is evidence for shear restructuring of 
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DLCA aggregates to larger fractal dimensions in colloids (Martin, Wilcoxon and Schaefer, 1990; Carpeneti, 1990; 
Jung, Amaland Raper, 1996; Selomulya, Bushell, Amal and Waite, 2002). There is, however, no evidence for ca. 
lOjltm , D ^ 1.8 aggregates in our flames necessary for restmcturing. A second scenario, that we consider more 

likely and have described previously (Sorensen, Kim, Fry, Shi and Chakrabarti, 2003; Kim, Sorensen and 
Chakrabarti, 2004; Fry, Chakrabarti, Kim and Sorensen, 2003), was inspired by our recent simulations (Fry, 
Chakrabarti, Kim & Sorensen, 2003) of 3d DLCA that show when the aggregates grow, they eventually fill all space 
and then jam together. The result is a percolated, D ~ 2.6 superaggregate of the smaller D — 1.8 , DLCA 
aggregates. Structure factors of our simulations, Fig. 4 is an example, show both the D — 2.6 and D ^ 1.8 regimes 
as seen in our flame data, e.g., Fig. 2a. We have shown that the crossover length scale between the D — 1.8 and 2.6 
regimes is consistent with soot volume fraction measurements (Sorensen, Kim, Fry and Chakrabarti 2003; Fry, 
Chakrabarti, Kim and Sorensen, 2003, Fry, Sintes, Chakrabarti and Sorensen, 2002). Figure 5 shows a diagram of 
our vision of a superaggregate. 

Figures 6a shows 7nsec exposure photographs of submillimeter soot in acetylene. These clusters are flat, 
quasi-2d (they are seen face on in Figs. 6a) similar to those seen previously in acetylene flames (Sorensen and 
Hageman, 2001) (stage 3). Following that previous analysis, cluster images were digitized into binary format (Figs. 
6b) then were analyzed to determine the number of dark pixels N p , which is proportional to cluster mass, and R g , the 

cluster radius of gyration for each cluster. If the clusters are fractal, then the well known relationship N p ~ R^ will 

hold. Figures 6c shows a log-log plot of N p vs. R g yielding a fractal dimension of D = 1.44 ± 0.05 . Thus flames 
yield flat aggregates with a fractal dimension of ca. 1.45 over length scales of ca. 10 to lOOjLtm . This morphology is 

stage 3 in Table 1. Similar results were found previously (Sorensen and Hageman, 2001). 

In previous work on laminar acetylene/air flames we concluded that this phase of soot occurs because the 
soot, which is confined to an annular flame front, becomes larger than the thickness of the flame front and thus is 
confined in an effectively two-dimensional (2d) space (Sorensen and Hageman, 2001) in which to aggregate. 
Simulations of 2d DLCA yield fractal aggregates with D=1.44 consistent with our measurements (Meakin, 1999). 
Thus we named this phase “two-dimensional soot.” Since this phase of soot is immediately larger than the second 
phase of soot described above, and since it has a different fractal dimension, 1.45 as opposed to 2.6, we propose that 
it is a third phase of soot consisting of 2d, DLCA aggregates of the superaggregates of phase two. We call these 
aggregates of superaggregates “supersuperaggregates .” They have three levels of structure, each level having its 
own fractal dimension, and two crossover length scales, one between each consecutive level, as outlined in Table 1, 
stages 1, 2, and 3. 

Flame spanning networks of soot are seen higher in the flames or at larger fuel flow rates . Example of a 
digitized image for acetylene flames is given in the lower left hand corner of Fig. 7. Once again we stress that as for 
the supersuperaggregates these networks are confined to the thin cylindrical annulus of the flame front. In previous 
work we argued that this network was an aerogel (Sorensen, Hageman, Rush, Huang and Oh, 1998), an argument 
that we will now support quantitatively. Analysis of this aerogel network structure was accomplished by Fourier 
transformation of the digitized images obtained with the photomicroscope. Analysis of Fourier transform data at 
low q is hampered by the finite size of the picture that is 1mm square. The raw data show “interference ripples” due 
to the hard edges of the picture box. To remove this box effect the uniform box was Fourier transformed, and the 
raw data were normalized with this to yield a corrected structure factor showing only the structure of the image, not 
the picture box. The resulting q -space function is a structure factor, S(q) vs. q, is proportional to a light scattering 
result, I(q), shown in Fig. 7. 

In each analysis at large q (> 1 000cm 1 ) a slope of -3 is found (not shown in the figures) We interpret this as 

Porod scattering, I (q) ~ q" (d4) for d=2, to imply that the image is uniformly black for lengths q" 1 < 10" 3 = 10|Lim . 

This value is consistent with our photographic resolution limit which we estimate to be 5pm , i.e., for length scale 
less than 5 to 10pm the image is uniformly black and two dimensional. This is a mild test of our analysis. 

For 200 < q < 700cm" 1 , which corresponds to length scales q' 1 of 14 to 50pm , the I(q) data roughly 
follow a power law of l( q) ~ q" 14 . This is a poor resolution observation of a remnant of the 2d, D — 1.45 
supersuperaggregates, stage 3. We propose that it is these supersuperaggregates that come together to form the gel 
network. For 30 < q < 200 cm" 1 , i.e., length scales of 50 to 300|Lim , data for all flames showed a I(q)~q" 2 1 power 

law. This power law can be understood by comparing to simulations of 2d DLCA aggregation we performed. 

These simulations were allowed to proceed until the system percolated to form a gel network. Figure 8 shows the 
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structure factor for a simulation with nearly 4 x 1 0 5 monomers allowed to aggregate via Brownian motion in a 2d 
box until only one large cluster formed. Real space analysis of this large cluster showed it had a fractal dimension 
of 1.9, consistent with a 2d, percolation cluster (Stauffer and Aharony, 1992). The entire box was Fourier 
transformed for direct comparison to the flame network data. The results in Fig. 14 are strikingly similar to the 
flame data including both the q" 14 and q" 21 regimes (and interference ripples due to the finite system size). The 

power law exponent (2.1) does not equal the fractal dimension, D-1.9, in this case because D is almost equal to the 
spatial dimension d. Note that when D=d, a Porod law with power -(d+1) results, which would be q" 3 for this 2d 
situation. Thus the q" 2 1 power law implies that the fractal dimension of the network, in both the flame soot and the 
simulation, is less than but very close to the spatial dimension of d=2. This strongly supports the contention that the 
gel-like network structure in our flames is the result of a 2d percolation of the D — 1.4 , 2d DLCA, 
supersuperaggregates . 
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